DESY 07-184 



Gravitino Dark Matter with Broken R-parity 

Alejandro Ibarra^ ^ 

Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany 



Abstract. Scenarios with gravitino dark matter face potential cosmological problems induced by 
the presence of the Next-to-Lightest Supersymmetric Particle (NLSP) at the time of Big Bang 
Nucleosynthesis (BBN). A very simple, albeit radical, solution to avoid all these problems consists 
on assuming that i?-parity is slightly violated, since in this case the NLSP decays well before 
the onset of BBN. Remarkably, even if the gravitino is no longer stable in this scenario, it still 
constitutes a very promising dark matter candidate. In this talk we review the motivations for this 
scenario and we present a model that naturally generates i?-parity breaking parameters of the right 
size to produce a consistent thermal history of the Universe, while satisfying at the same time all 
the laboratory and cosmological bounds on i?-parity breaking. We also discuss possible signatures 
of this scenario at gamma ray observatories and at colliders. 
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1 Motivation 

The gravitino, when it is the lightest supersymmetric 
particle, constitutes a very promising candidate for the 
dark matter of the Universe. The interactions of the 
gravitino are completely fixed by the symmetries, and 
the thermal relic density is calculable in terms of very 
few parameters, the result being [1] 

(1) 

while the relic density inferred by WMAP for the ylCDM 
model is fl^^i^^? — 0.1 [2]- In this formula, Tji is the re- 
heating temperature of the Universe, is the grav- 
itino mass and m~ is the gluino mass. It is indeed 
g ° ^ 

remarkable that the correct relic density can be ob- 
tained for typical supersymmetric parameters, 7713/2 ^ 
100 GeV, m~ ~ 1 TeV, and a high reheating temper- 
ature, Tr ~ 10^° GeV, as required by baryogenesis 
through the mechanism of thermal leptogenesis 3 . 

If thermal leptogenesis is the correct mechanism 
to generate the observed baryon asymmetry, then the 
NLSP can pose serious cosmological difficulties. If one 
assumes exact i?-parity conservation, as is commonly 
done in the literature, then the NLSP can only decay 
into gravitinos and Standard Model particles with a 
decay rate strongly suppressed by the Planck mass, 
yielding a lifetime 
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The high reheat temperature for the Universe re- 
quired by thermal leptogenesis, Tr > 10^ GeV [4], and 
the requirement of correct relic abundance, Eq. ([T]), im- 
ply m3/2 ^5 GeV for a gluino mass of m~ — 50GGcV. 
Therefore, the lifetime of the NLSP is typically of sev- 
eral days and it is present at the time of BBN, po- 
tentially jeopardizing the successful predictions of the 
standard scenario. Indeed, this is the case for most 
supersymmetric scenarios with gravitino dark matter. 
Namely, when the NLSP is a neutralino, its late decays 
into hadrons can dissociate the primordial elements ^ , 
and when the NLSP is a stau, its presence during BBN 
catalyzes the production of ^Li, yielding a primordial 
abundance in stark conflict with observations [6j. 

One simple, albeit radical, solution to the problems 
induced by the NLSP is to assume that i?-parity is not 
exactly conserved [7 . In fact, there is no deep theoret- 
ical reason why i?-parity should be exactly conserved, 
and although present experiments indicate that R- 
parity is approximately conserved, a slight violation of 
it cannot be precluded. Without imposing any ad hoc 
discrete symmetry, the superpotential of the Minimal 
Supersymmetric Standard Model (MSSM) reads [S] 

W = Wrp + ^\,,kUL,el + A^;, UQ.dl 

+ \^'Ijkutd''jdl + fi,L,Hu , (3) 

where Wnp is the familiar superpotential with i?-parity 
conserved. 

When i?-parity is broken, new decay channels are 
open for the NLSP, apart from the strongly suppressed 
gravitational decay into gravitinos and Standard Model 
particles. Namely, when the NLSP is mainly a right- 
handed stau, it can decay tr ^Vr, through the 
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coupling A323, with lifetime 



lO^s 



A: 



323 



10 



-14 



100 GeV 
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Thus, even with a tiny amount of i?-parity violation, 
A > 10""'^^, the stau NLSP will decay before it can 
have any significant impact on Big Bang Nucleosyn- 
thesis. A similar argument can be applied for the case 
of neutralino NLSP with analogous conclusions. 

On the other hand, there are very stringent labora- 
tory and cosmological constraints on the size of the R- 
parity breaking parameters. The strongest bound for 
the heaviest generation stems from the requirement of 
successful baryogenesis. Since A, A' and A" violate lep- 
ton and baryon number, if they are on thermal equi- 
librium at the same time as the sphaleron processes, a 
preexisting baryon asymmetry will be erased. There- 
fore, the requirement that these couplings are out of 
thermal equilibrium before the electroweak phase tran- 
sition translates into the bound [S] 



A, A', A" < IO- 
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This bound is a sufficient but not a necessary condi- 
tion, and could be relaxed for some flavour structures. 

Interestingly, if the size of the i?-parity violating 
couplings is in this range, 10"^'' ^ A, A' < 10~^, the 
gravitino still constitutes a viable dark matter can- 
didate [To]. The decay rate is doubly suppressed by 
the inverse Planck mass and by the i?-parity breaking 
couplings, yielding a lifetime 
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which is several orders of magnitude longer than the 
age of the Universe. 

In the next sections we will present a model that 
can naturally accommodate our favoured range of im- 
parity breaking parameters, 10"^'' ^ A, A' < lO"'', and 
we will discuss the possible experimental signatures of 
this scenario, both at gamma ray observatories and at 
colliders. 



2 A model for small (and peculiar) 
i?-parity violation 

We will use for convenience 5*0(10) notation, although 
our model is not necessarily embedded into a Grand 
Unified group. Then, quarks and leptons will be de- 
noted by 16i and the Higgses by IOh. In order to give 
Majorana masses to neutrinos B — L has to be broken, 
either by the vacuum expectation value of a 16, 16 
(with B - L = ±1) or of a 126 (with B - L = 2). To 
have only particles with small representations we will 
break B — L with a 16 and a 16. Thus, i?-parity will 
be necessarily violated after the spontaneous breaking 
of the B — L symmetry at the scale vb-l- 

With this matter content there are two kind of 
terms in the superpotential allowed by the 50(10) 



symmetry. The terms 16il6jlOH and jgjl6il6jl6 16 

are "good terms" , since they produce Dirac masses and 
a Majorana mass term for the right-handed neutrinos, 
respectively. On the contrary, the terms 16il610H and 
jg^l6il6jl6kl6 are "bad terms"; the first one pro- 
duces after B — L breaking the i?-parity violating bi- 
linear term vb^lLHu that in turn produces too large 
neutrino masses, while the second one produces the im- 



parity violating trilinear terms 



-QLd" 



Mp Mp 

that induce too rapid proton decay. It is interesting to 
note that both "bad terms" are generated when the 
16 representation acquires a large vacuum expectation 
value. Since the existence of the 16 representation in 
the spectrum is unavoidable, the only way to suppress 
these couplings is by means of additional symmetries. 

To this end, we will impose in our model a U{1)r 
symmetry and the following assignment of charges: 





16, IOh 16 16 


1 


R 


10 0-2 


-1 



where the singlet 1 has been introduced in order to 
break the i?-symmetry. With this assignment of charges, 
holomorphicity guarantees that after B — L breaking 
no i?-parity violating term will be generated from the 
superpotential at any order in perturbation theory. 
Nevertheless, the Kahler potential is not holomorphic 
and could be a source of i?-parity violation. Indeed, 

terms such as 1 16^16^10//, 1^1616,10// can appear 
in the Kahler potential, and will eventually produce 
bilinear i?-parity violation. 

With these elements in mind, it is relatively sim- 
ple to construct a model that produces small lepton 
number violation and tiny baryon number violation. 
In Standard Model notation, the previous model is: 





Q, u^, ii„, Hd N 
d", L, ly" 


^ Z X 


B-L 


±1/3, ±1 1 -1 





R 


1 0-2 
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where the spectator fields <P, Z and X have been intro- 
duced in order to break B — L (and i?-parity), to break 
supersymmetry, and to ensure that {N) = (N'^) = 
vb~l, respectively. After these singlets acquire a vac- 
uum expectation value, the effective superpotential reads 
W =i Wd + Wy. + Wo . Here, Wb is the i?-parity 

conserving superpotential of the MSSM and W^^ is 
the part of the superpotential involving right-handed 
neutrinos. In this model the heaviest right-handed neu- 
trino mass is given by M3 ~ Vb-l/Mp, which is nat- 
urally large, leading through the see-saw mechanism 
to small neutrino masses. Finally, Wj^ is the i?-parity 

breaking superpotential given in Eq. ([3]) . We will choose 
to work in the Li — Hd basis where fii — and thus 
all the i?-parity violation is encoded in the trilinear 
terms. In this basis, 
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The lepton number violating couplings are sup- 
pressed compared to the charged lepton and down- 
type quark Yukawa couplings, /i^ and h''- respectively, 
by v'^b-lI-^p ^'^'i by some coefficients C = 1.0. ..0.01 
that depend on the flavour structure of the Kahler po- 
tential. On the other hand, the baryon number violat- 
ing coupling only arises after supersymmetry breaking 
and is suppressed by higher powers of the Planck mass. 
To estimate the scale oi B — L breaking and the size 
of the coefficients C, we will use the flavour model 
proposed in [11]. For this particular model, we obtain 
A3ij,A3jj ~ 10~®, within the desired range 10"^"* < 
A, A' < 10^^, and A" ~ 10~^^, yielding negligible rates 
for proton decay [7]. 



3 Signatures at gamma-ray observatories 

When i?-parity is broken, gravitino decays could be 
happening at fast enough rates to make its decay prod- 
ucts detectable in future experiments. In particular, 
the photon flux produced by the gravitino decay could 
be detected as an extragalactic diffuse gamma-ray flux 
with a characteristic spectrum, that is originated from 
the decay of gravitinos at cosmological distances and 
from the decay of gravitinos in the Milky Way halo. 

In the case that the gravitino is lighter than the 
W boson, the gravitino decays mainly into a photon 
and a neutrino, producing a photon spectrum consist- 
ing in a monochromatic line. If gravitinos constitute 
the dominant component of dark matter, the decay of 
gravitinos at cosmological distances will be detected at 
the Earth as a perfectly isotropic diffuse gamma ray 
background, with an energy spectrum corresponding 
to a red-shifted monochromatic line: 



dE 



1 



3 

flu 



-1/2 



2/^/2 0(l-y) , (8) 



where y ~ lEjm-j^ii and 

^Zl2Pc _ 10-6 GcV / T3/2 ^ -1 



1/2 



cm^ str s V 10^'^ s 



(9) 



Here, ^23/2^^ = 0.1, = 1.05 h'^ x IQ-^GeVcm-^, 
= 0.25, Ha ^ h 100 km s'^ Mpc"^ with h = 
0.73 [H] and T3/2 is given by Eq. 

In addition to the contribution to the photon flux 
from the decay of gravitinos at cosmological distances, 
one also expects a monochromatic line stemming from 
the decay of gravitinos in the Milky Way halo. This 
contribution reads: 



I dj halo 

dE 



2E 



where 



1 



Phalo{l)dl 



(10) 



(11) 



13/2 Jl.o.s. 

Since the Earth is not located at the center of the 
Galaxy, the gamma-ray flux from the decay of grav- 
itinos in the halo is slightly anisotropic. For typical 



halo models, we find that the halo component domi- 
nates over the cosmological one, giving rise to a slightly 
anisotropic gamma ray flux with an energy spectrum 
dominated by a monochromatic line. 

There are very stringent limits on the photon flux 
produced by decaying dark matter from the observa- 
tions of the Energetic Gamma Ray Experiment Tele- 
scope (EGRET) aboard the Compton Gamma Ray 
Observatory. Assuming that one understands the galac- 
tic foreground, one can extract from the EGRET data 
the extragalactic diffuse component. The first analysis 
by Sreekumar et al. [13j gave a roughly isotropic ex- 
tragalactic flux with an energy spectrum described by 
the power law 



E^^ = 1.37 X 10-6 
dE 



1 GeV/ 



GcV 

cm^ str s 



(12) 



in the energy range 50 MeV-10 GeV. The improved 
analysis of the galactic foreground by Strong et al. [M] , 
optimized in order to reproduce the galactic emission, 
shows a power law behavior between 50 MeV-2 GeV, 
but a clear excess between 2-10 GeV, roughly the 
same energy range where one would expect a signal 
from gravitino decay. In view of all the systematic 
uncertainties involved in the extraction of the signal 
from the galactic foreground, we find it premature 
to attribute this excess to the gravitino decay. We 
find nevertheless this coincidence as interesting and 
deserving further attention. The upcoming satellite- 
based gamma ray experiments GLAST and AMS-02 
will measure the energy spectrum with unprecedented 
accuracy, providing very valuable information for the 
scenario of decaying gravitino dark matter. 

We show in Fig[l]the expected signal for a decay- 
ing gravitino with a mass of 10 GeV and a lifetime of 
10^^ s [H]. To compare our results with the EGRET 
data [11], also shown in the flgure, we have averaged 
the halo signal over the whole sky excluding a band of 
±10° around the Galactic disk, and we have used an 
energy resolution of 15%, as quoted by the EGRET 
collaboration in this energy range. The photon spec- 
trum is dominated by the sharp line coming from our 
local halo, while the red-shifted extragalactic compo- 
nent is somewhat fainter. 

If the gravitino is heavier than the W or Z bosons, 
new decay channels are open. In this case, the energy 
spectrum consists of a continuous component, stem- 
ming from the fragmentation of the gauge bosons, and 
a relatively intense gamma ray line [16j . 



4 Signatures at colliders 

In the scenario proposed in this work, the lifetime of 
the NLSP is around one nanosecond or longer, giving 
rise to very distinctive signatures at colliders. When 
the NLSP is the lightest stau, that we assume mainly 
right-handed, the main decay channel is tr — > r z^^ , /i ^-7- 
The corresponding decay length is 



ct]^p ~ 30 cm ( 



200GeV 
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(13) 
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Fig. 1. Energy spectrum of extragalactic (long-dashed 
green) and halo signal (dotted magenta) compared to the 
EGRET data, for rg/a = lO^'^ s and mg/a = 10 GeV. 
The data points are the EGRET extragalactic background 
extracted by Strong et al. in [14], while the short-dashed 
(blue) line shows the powerlaw fit from Sreekumar et al, 
Eq. (|12|l . Finally, the solid (thick red) line shows the sum 
of these contributions with a powerlaw background (thin 
red line), which has been obtained fitting the low energy 
EGRET points. 



a small amount of i?-parity violation, 10~^^ < A, A' < 
10"^, yields a consistent thermal history of the Uni- 
verse, in the sense that it allows baryogenesis through 
thermal leptogenesis, provides a viable candidate for 
cold dark matter, and does not spoil the successful 
predictions of the Standard Big Bang Nucleosynthesis 
scenario. 

We have presented a model that links i?-parity 
breaking to B — L breaking, and generates i?-parity 
breaking parameters in the right range to accommo- 
date the abovementioned scenario. The relic gravitino 
decays into photons produce a diffuse halo and extra- 
galactic gamma ray flux, that might have been ob- 
served already by EGRET. Future experiments, such 
as GLAST or AMS-02 will provide unique opportu- 
nities to test the scenario of decaying gravitino dark 
matter. Finally, we have discussed the striking signa- 
tures that this scenario might produce at future col- 
liders, consisting in a vertex of the NLSP significantly 
displaced from the beam axis. 

Acknowledgments. I would like to thank my col- 
laborators G. Bertone, W. Buchmiiller, L. Covi, K. 
Hamaguchi, D. Tran and T. T. Yanagida for a very 
pleasant and fruitful collaboration. 



Besides, the small left-handed component of the stau 
mass eigenstate can trigger a decay into jets through 
Tx — > bH, provided the process is kinematically open. 
The hadronic decays are enhanced compared to the 
leptonic decays by the larger bottom Yukawa coupling 
and by the colour factor, but are usually suppressed 
by the small left-right mixing. The decay length in this 
channel reads 

(14) 

where 9^ is the mixing angle of the staus. Therefore, 
if the decay occurs inside the detector, the signature 
would consist on a heavily ionizing charged track fol- 
lowed by a muon track, by one jet or by three jets. 

On the other hand, when the NLSP is the light- 
est neutralino, it decays through Xi ~^ t^W^ if this 
decay channel is kinematically accessible [17j . The cor- 
responding decay length can be approximated by 

producing a significantly displaced vertex followed by 
jets, that could be observed provided the decay occurs 
inside the detector. If this decay channel is kinemat- 
ically closed, the neutralino typically decays outside 
the detector yielding identical collider signatures to 
the case with i?-parity conserved. 



5 Conclusions 

We have argued that a supersymmetric scenario with 
gravitino LSP with a mass in the range 5-100 GeV and 
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